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Abstract 



The strong Scott conjecture about the electron density at a distance 
1/Z from an atomic nucleus of charge Z and its generalization for molecules 
are proved. The density, suitably scaled, converges to an explicit limiting 
density as Z — > oo. Both a weak convergence and a pointwise convergence 
on spheres holds. 

1 Introduction 

A great deal is known about the ground states of large atoms in the framework 
of the non-relativistic Schrodinger equation, with fixed (i.e., infinity massive) 
nuclei. The leading term, in powers of the nuclear charge Z, is given exactly 
by Thomas-Fermi theory, as was proved by Lieb and Simon fll2| |; see JT]] | for 
a review. This leading term in the energy is proportional to Z 7 / 3 , with the 
proportionality constant depending on the ratio of N/Z, which is assumed to 
be held fixed as Z — ► oo. Here, N is the electron number. Neutrality, i.e., 
N = Z, is not required, even though it is the case of primary physical interest. 
The characteristic length scale for the electron density (in the sense that all 
the electrons can be found on this scale in the limit Z — ► oo) is Z~ 1//3 . The 
fact that the true quantum- mechanical electron density, p ( i, converges (after 
suitable scaling) to the Thomas- Fermi density, p TF , as Z — > oo with N/Z fixed 
was proved in |12| . The chemical radius, which is another length altogether, is 
believed, but not proved, to be order Z° as Z — > oo. 

The first correction to the Z 7 / 3 law is not, as was formerly supposed, the 
^5/3 corrections arising from exchange and correlation effects and kinetic en- 
ergy corrections on the Z -1 / 3 scale. Instead it is Z 2 = Z 6 / 3 and arises from 
extreme quantum mechanical effects on the innermost electrons, which are at 
a distance Z" 1 from the nucleus. Among these the most important are the K 
shell electrons. It was Scott |l9| who pointed this out and he also gave a formula 
for the correction term to the energy, 
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where q is the number of spin states per electron (of course q = 2 in nature). 
It is noteworthy that E Scott does not depend on the fixed ratio N/Z, provided 
N/Z ^ 0. This fact agrees with the idea that E Scott arises from the innermost 
electrons, whose energies, in leading order, are independent of the presence of 
the electrons that are further from the nucleus. The truth of (1) (i.e., the 
statement that the energy is E TF + E Scott + o(Z 2 ) for fixed N/Z ^ 0) was 
proved in ||| (upper and lower bound) and by Hughes || (lower bound) 
in the atomic case and by Bach § in the ionic case. (For different proofs and 
extensions of this result see , Fefferman and Seco || 0, Q 0, R|, and Ivrii 
and Sigal |l(J.) 

The "strong Scott conjecture", which we prove here, was made later in jyj, 
Equation (5.34). It concerns the electron density pd at distances of order Z^ 1 
from the nucleus and states that in limit Z — ► oo a suitably scaled pd converges 
to the sum of the squares of all the hydrogenic bound states. This function 
is defined in Section 2 below and is extensively analyzed in (Previously, 
an upper bound for pd at the origin of the correct form, namely 0(Z 3 ), was 
derived in [^0[ |l[.) We prove this convergence in several senses, one of which 
is a "pointwise" convergence on spheres. In fact we go further and show that 
the individual angular momentum densities converge to the hydrogenic values, 
thereby giving a somewhat more refined picture of the ground state. 

Our main proof strategy is the usual one. We add e times a one-body test 
potential to our Hamiltonian and then differentiate the ground state energy with 
respect to e at e = in order to find pd- To obtain pointwise convergence the 
test potential is a radial delta-function. To control the energy we rely, in part, 
on the results and methods in |22, 23 



In the following we shall state and prove our theorems for the neutral case 
N = Z. We do so to avoid the notational complexity and additional discussion 
required for N/Z ^ 1. It is straightforward, however, to generalize our results 
to N/Z ^ 1. 

In the next section precise definitions, as well as our main theorems are 
given. Section 3 contains two lemmata about the difference in energies with and 
without the test potential. Since there are infinitely many hydrogenic bound 
states, we need these estimates in order to be able to show that the sum of 
the derivatives (with respect to e) equals the derivative of the sum. The strong 
Scott conjecture for atoms is proved in Section 4. Section 5 contains the obvious 
extension to molecules. The Appendix contains a few needed facts about ground 
state energies. 

2 Definitions and Main Results 

The Hamiltonian of an atom of N electrons with q spin states each and a fixed 
nucleus of charge Z located at the origin is given by 



N 

Hn.z = 



in units in which h 2 /2m — 1 and |e| = 1. It is self-adjoint in the Hubert space 

JV 

fjN '■= A (£ 2 (H 3 ) C 9 ), i.e., antisymmetric functions of space and spin. A 
v=x 
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general ground state density matrix, denoted by d, can be written as 



M 



d = '^2w v \ip I/ }{ip I/ \, (3) 

where the ijj v constitute an orthonormal basis for the ground state eigenspace 
and where the w v are nonnegative weights such that X^^Li w v = 1. It is well 
known that the ground state can be degenerate, e.g., it is for the carbon atom. 
The corresponding one-electron density is the diagonal part of the one-electron 
density matrix and is, by definition, 

M q . 

p d (v)=Ny^w v y~] / \tp„(x,V2,---,x N ;<7 1 ,...,a N )\' 2 dx 2 ---dv N . 

1 , JR3<«-1) fA \ 

v=l ai,...,a N =l " [Jij 

The density p^d of angular momentum I electrons at radius r from the nucleus 
is given in terms of the normalized spherical harmonics Y/ iTO (w). 

I M 

Pi 



/ F; m (w)V'i/(rw,t2, ...,tjv;o-i, ...,a N )dio 
Js 2 



2 

dx 2 .--dx N (5) 



where we write r = ru; and denotes the usual unnormalized surface measure 
on the two dimensional sphere S 2 with (47r)^" 1 J §2 cL; = J §2 |F; jm | 2 da) = 1. 

Throughout the paper we will write Lp T z F (r) for the Thomas- Fermi potential 
of electron number N — Z and nuclear charge Z, i.e., 

v¥(r) = Z/r-\t\- 1 *pY, 
where p T z F is the nonnegative minimizer of the Thomas-Fermi functional 

(f (6tt 2 /<z) 2/ V /3 (r) - ^ P (x))dx + D(p, p) 

under the condition J p = N = Z and with 

D(p,p):= l -l ^Ml dxds . 



Both ipz and p T z F are spherically symmetric, i.e., they depend only on r = 
|t|. There is a scaling relation (p^ F (r) = Z i ^ip TF '(Z 1 / 3 r), where (pJ F is the 
Thomas- Fermi potential for Z = 1 and "electron number" equal to 1. Similarly, 
Pz F ( r ) = Z 2 p[ F {Z 1 / 3 r). This scaling shows that the "natural" length in an 
atom is Z -1 / 3 . Note that the Thomas- Fermi functional has a unique minimizer. 

The Scott conjecture, on the other hand, concerns the length scale Z^ 1 , 
where we expect the density to be of order Z 3 instead of Z 2 . In terms of the 
"true" density defined in (4), we now define 

Pz (x) := Z- 3 p d (x/Z). (6) 
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Likewise, we define the angular momentum density 



p LZ (r) := Z- 3 p lA (r/Z). (7) 

To formulate the strong Scott conjecture we consider the angular momentum 
I states of a hydrogen atom (Z — 1) with radial Hamiltonian 

»f * + !2+i>_i (8, 
ar z r z r 

with normalized eigenfunctions ipf 1 (that vanish at and oo) corresponding 
to negative eigenvalues ef n . (The superscript H denotes "Hydrogen" and dis- 
tinguishes hi from other radial Hamiltonians to be considered later.) The 
normalization is j Q \if} n ^(r)\ 2 dr — 1. We define the corresponding density in 
the channel I to be 

oo 

pf (r) := q(2l + K(r)| 2 /(47rr 2 ); (9) 

n=0 

the total density is then 

oo 

p H (r)=J2P?(.r). (10) 
i=o 

Although we shall not be interested in detailed properties of p , we note 
the following proved in The sum over I and n defining p H (r) is pointwise 
convergent for all r. It is monotone decreasing and it decays asymptotically for 
large r as l/(6-7r 2 r 3 / 2 ). This large r asymptotics meshes nicely with the small r 
behavior of ^pJ F (r). 

Note: In [|| the operator h H is defined using atomic units H 2 /m = 1, i.e., 
with \{~d 2 /dr 2 + 1(1 + l)/r 2 ) instead of -d 2 /dr 2 + 1(1 + l)/r 2 . Note also that 
we have included the factor q in the definition of p which was not done in 
Thus some care is needed in comparing formulae there with formulae here. 

Various notions for the convergence of the rescaled density are possible. Our 
precise statements are Theorems [l] and ^ below and Theorem ^ in Section g. 

Theorem 1 (Convergence in angular momentum channels). Fix the 

angular momentum Iq- 



1. For positive r 



(pointwise convergence). 



hm p^ z (r) = p?(r) (11) 

Z— *oo 



2. Let V be an integrable functions on the positive real line. Then we have 
the weak convergence 

lim / rV(r)p kuZ (r)dr = / rV(r)p"(r)dr. (12) 
z^°°Jo Jo 
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Theorem 2 (Convergence of the total density). 

1. Let W be a bounded (not necessarily constant) function on the unit sphere. 
Then, as Z — > 00, 

/ W{ui)p z (ru))dw -> p H {r) [ W(iv)duj (13) 
Js 2 Js 2 

(pointwise convergence of spherical averages). 

2. Let V be a locally bounded, integrable function on R 3 . Then, as Z — ► 00, 

/ \t\V(t) Pz (x)dt^ I \x\V(t)p H {\x\)dt. (14) 

Remarks: 

1. It is not really necessary to take a sequence of ground state density matri- 
ces. We could take just a sequence of states, dN,z, that is an approximate 
ground state in the sense that 

tr (HN,zdN,z) — Ejy^z „ 
Z 2 

as Z — > 00. Here En,z is the bottom of the spectrum of Hn,z- It might 
not be an eigenvalue, and it certainly will not be one if N/Z is larger than 
2. 

2. It is important to note that W and V in (|l^) and (|14 ) need not be spher- 
ically symmetric. It might appear that only the spherical averages of W 
and V are relevant, but this would miss the point. Theorem^ says, that in 
the limit Z — > 00 there is no way to construct a ground state or approxi- 
mate ground state that is not spherically symmetric on a length scale Z . 
For example, in the case of carbon there are ground states that are not 
spherically symmetric and for which replacing W by its spherical average 
changes the left side of jl^). 

3. A word about pointwise convergence. The one-body density matrix j(t, t'), 
which is defined as in but with \ipi/(-- -)| 2 replaced by 

V>„(t, r 2 , ...,tjv;01,.. ■ ,<7jv)^(t',r 2 , ...,tjv;0l,.. .,o"jy), 

is easily seen to be in the Sobolev space i? 1 (R 3 ) when 7 is considered 
as a function of each variable separately. The trace theorem in Sobolev 
spaces then implies that the function of u> on the sphere § 2 , j(rui, r'ui') is 
in L q (S> 2 ) for all q < 4. Thus, the integrals in (4) are well defined and 
Pd(ruj) — "f{ru>,ruj) is in i 2 (S 2 ). It is also easy to see that y/ pi,z(r) is 
in i? 1 (0,oo) and hence it is a continuous function of r. Since pd(x) is in 
i 2 (S 2 ) the integrals in (|ll|) and (g) are well defined when W S i 2 (S 2 ). 
— If p and 7 belong to a ground state of H^^z with N < Z then they are 
even continuous functions in all variables. This follows from the regularity 
theorem of Kato and Simon (Reed and Simon |l8| , Theorem XIII. 39) and 
the uniform exponential decay of ground state eigenfunctions. This decay 
is implied as follows. By Zhislin's theorem the atomic Hamiltonian has 
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infinitely many eigenvalues below the essential spectrum and the ground 
state eigenspace has finite dimension. This implies that the ground state 
energy is always a discrete eigenvalue which, in turn, implies exponential 
decay of the ground state eigenfunctions according to Theorem XIII. 42 of 




3 Eigenvalue Differences of Schrodinger Opera- 
tors Perturbed on the Scale 1/Z 

It is well known, and will be seen more explicitly in Section 4, that the eigenval- 
ues of H]y t z can be controlled to within an accuracy of o(Z 2 ) by considering a 
one-body Schrodinger operator with the spherically symmetric potential given 
by Thomas-Fermi theory. In the angular momentum I channel, this is 

tff = + ^ - *S'(r). (15) 

(We suppress the dependence on N in hJ F since N = Z.) Closely related to 
hJ F is the unscreened hydrogenic Hamiltonian 



In this section we want to study how the spectra of these operators are 
shifted by the addition of a perturbing potential of the form 

eU z (r) = eZ 2 U(Zr) 

where e is a small parameter and where U is some fixed function. In particular, 
U will be a radial delta function, U(r) — S(r — a) for some a > 0. 
Both cases, h TF and h H , will be considered together and we write 

hi, e ,z = +V ilZ (r) - eU z (r), (17) 

in which V lt z = -Z/r + l(l + l)/r 2 or V uz = -(p% F (r)+l(l + l)/r 2 . 

Our first lemma estimates the difference in the spectra of hifl,z and hi >e> z 
by the difference in the trace (tr) of the negative parts (hi :<L: z)- or (hifi,z)- 
(i.e., the sums of the negative eigenvalues). This lemma will later on allow us 
to interchange the limits Z — ► oo and e — * with the I summation. 

Lemma 1 Set U(r) = S(r — a), Uz(r) = Z 2 U{Zr) and assume \e\ < n/(16a). 
Then 

| trCfcy,,*)- - tr(M**)_ I < ^ i + 1 9 ) °f 2Z + 1) - 
Proof: By the minimax principle we have for e > 

< s e< i >z := ti(hi fi ,z)- - ta(hi, E ,z)- < etr(U z di ieiZ )- (18) 
Inserting the identity twice in the right side of we have 

s e ,i,z < etr (ABCB* A*) < e\\ AfJjBf^ trC (19) 
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with 



B 

C 



= (hi, e , z + d,z)- 1/2 {H ,i + ci,z) 1/2 , 

= (H ,i + c^ z r 1/2 U z (H Q j + c ltZ y 1/2 > 0, 



where c^z is any positive number bigger than | inf tr(/i; jCj z)|, where a(h) denotes 
the spectrum of h. We also define H i := —d 2 /dr 2 + 1(1 + l)/r 2 to be the free 
operator in the angular momentum channel I. Since <p% F (r) < Z/r and since 
inf a(H ,i - Z/r) = -Z 2 /[A(l + l) 2 ] we can take c/, z := Z 2 /(l + l) 2 provided e 
is not too large. 

We now estimate these norms individually: 

Because c^z is bigger than the modulus of the lowest spectral point of ft-z, e ,z 
and di^_^z is the projection onto the negative spectral subspace of /i;, e ,z we have 

Plloo < V^- (2°) 

For B we get 

\\B4>\\ 2 = \\(hi, e , z + %zT 1/2 {H ,i + Q,z) 1/2 ^|| 2 

= (</>, (Ho, + Q, z ) 1/2 (Ju, e ,z + c hz )-\H 0il + c hZ ) 1/2 4>) 

-^db* (21) 

with 

Wi, e , z := (H ,i + Q, z )- 1/2 &z + eUz)(H Q ,i + Q,z)" 1/2 . 
We will then have 

\\B\\<V2 (22) 

if we can show that Wi >e> z is bounded above by |. To this end we note that 
Hq,i + c^z is invertible, so that we can write any normalized G L 2 (R + ) as 
4> := (Ho,! + CLz) 1/2 ip/\\(H ,i + ctz) 1/2 ip\\ with V in the domain of H ,i. Thus, 
we have to show that 

(^ Wj, e ,z$ = (V, faz + eU z )ip)/(i>, (H ,i + c ltZ m < \, 
which is equivalent to 

^(H ,i + ci, z ) - <Pz - eU z > 0. (23) 
Since (p^ F (r) < Z/r and since 

,-a/Z ,-a/Z 

(tp, U z ip) = Z\^(a/Z)\ 2 <Z \4>\ 2 '(r)dr < 2Z3? / tp(r)tp' (r)dr 



a/Z 



1/2 



<2Z\W\\A / ^(r) 2 dr\ (24) 
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we have that (tp,Uztp) < (4a / 7r) 1 1 -t/>' Hi- (Here we use the inequality, J Q L -tp' 2 > 
(ir/2L) 2 j£ ip 2 , when ^(0) = 0.) Thus, @ is implied by 

- 4 i7^ + «.* + ^ v w * °- (25) 

The sum of the first two terms in ( p5| ) vanishes because of our choice of ci t z and 
last term in ( p5| ) is zero when e < 7r/(16a). This is true by hypothesis, and the 
bound on the norm of B is proved. 

Finally the trace of C is computed easily, since it is of rank one. Since the 
kernel {Hqj + ci t z) [r, r'), is a positive, continuous function in both variables, 
tr C = Z(H + ci,z)~ 1 {a/Z, a/Z). A well known calculation yields 

{Ho,l + ci l z)~ 1 {r,r') = <frK l+ i(y/ciJ r>)I l+ 1 {J%E r<)V7, 



where r> = max{r, r'} and r< = min{r, r'}. Thus 

trC = aK l+ i(y/c[^—)I l+ %(y/c[^—). 

The modified Bessel functions I l+ i and K l+ i are both positive and the following 
uniform asymptotic expansions hold. (See Olver for a proof of the estimates 
of the remainder terms, fl5|| , p. 6 for the remainder in the form used here, |l7jj , 
Chapter 10, Paragraph 7 for a review; see also Olver fl6| , section 9.7.) 

K n (nx) = J e" n « [l + e ,2 (n,t)] (26) 



where 



t e n£ - 

In{nx) = — - —-[l + e 0;1 (n,t)], (27) 

" zttti 1 — eo,i( n , 0) 



£ : = T - 7T lo S ■ 



and 



t 2 1 -t 
t := (l + x 2 )~5 



|eo,i(M)l < U ° 

n — uq 

n 



|eo,2(n,t)| < 

rt — no 

with n := -^m + < 1/6- Thus, 

6v5 -L* 

K n (nx)l n (nx) < r - — — < — , 

2n(l + x 2 ) 2 (n - 2n )(n - n ) in 

where the last inequality holds for n > \. Thus 

« c< -lmTTy (28) 



Putting (§§), and © together yields 

9q z a 9aZ 2 
s c ,lz < e < e- 



(2Z + 1) " (Z + l) 2 (2^ + l) 

which is more than the desired result for e > 0. 

If e is negative we have, again by the minimax principle, 

> s et i, z := tr(hi^z)- - tr(h he , z )- > etr(U z di A z). (29) 

Similar to the previous analysis, we have 

Se,l,z > etx(DCD*) > e||D||^trC (30) 

with 

D := di fi .z(hi fl ^ z + %z)' s ■ 

As for A above, we have 

||-D||oo < Vo^i". 
Putting this together with (^8|) gives 

^ 9c LZ a 9aZ 2 



'2(21+1) " 2(Z + 1) 2 (2Z + 1) 

which is even better than the desired result for negative e. | 

The next result will later on allow us to interchange the limits Z —> oo and 
e — ► with the n summation for fixed I. 

Lemma 2 Set U(r) = S(r — a) and assume \e\ < 7r/(4a), a > 0. Let 

d 2 1(1 + 1) 1 

/l ^ := -TT + ^- i eC/ ( r ) ( 31 ) 

ar r r 



\e n ,i,o - e ra ,;, e | < ; n2 777 (32) 



with form domain Hq(0,oo). Let e nt i >e denote the n-th eigenvalue of hi ie . Then 

1 ea 
(n + Z) 2 7T — 4ea 

Proof: For any -0 in Hq(0, oo) we have 

hMa)| 2 <2-aW\\l 

7T 



as proved in (24) of Lemma 1. Thus, for e > 0, 



4ea. 

hl,e > (1 ) 

7T 



<ir 2 r 2 (1 



(33) 



This implies 

, 4ea 
e„,i, e > (1 )e n ,i,o 

7T 



where e^^o is the n-th eigenvalue of [ ] in (p3|), i.e., where the potential r 1 is 
replaced by (1 — 4ea/7r)~ 1 r~ 1 . Thus, 

1 / „ .„ 4ea, A 1 ea 

< e„,i, - e n ,i,e < ., ■ -1 + (1 ) _1 



4(n + 1) 2 \ 7T / (n + Z) 2 7T - 4ea 
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which proves the claim when < e < ir/ (4a). 
If e is negative we have 



4ea, 

h, t < (1 ) 



dr 2 



1(1 + 1) 



(l-^> 



which again proves the claim (by the same argument) when > e > — ?r/(4a). 



4 Proof of the Strong Scott Conjecture 

We are now able to give the proof of our theorems. We begin with the proof of 
Theorem [j] and begin with the first statement: 

1. The proof of the convergence of the spherical averages: Set U(r) := S(r—a) 
for a > and U z (r) := Z 2 U{Zr) = ZS{r - §). Fix l and let 

N 

H' N>Z := H N . Z -eJ^UziruMlo) (34) 

v=l 

where II(Zo) denotes the projection onto angular momentum Iq. Wc define \{Z) 
- which does not depend on e - by 

X(Z) := a pi (a) = — 2 . (35) 

Let us define e i ^ l z and e n ,la,Z\ n = 1» 2, £ € 1, to be the negative eigen- 
values of the operators 

H l,e,z ■= + — j-eU z Si, k (36) 

d 2 1(1 + 1) TF 
Ht, e , z := -— + ±- s -L-^-eU z 5 l>h (37) 

with zero Dirichlet boundary on (0, oo). To obtain an upper bound for \(Z) we 
pick e positive and estimate as follows: by (^) we have the upper bound 

tr(H NtZ d) 

L — X oo 

< ^q(2l+l)^e* tlfitZ +^q(2l+l)^e n<lfi , z -D{ P TF,pTF)+covstZ& 

1=0 n l=L n (38) 

where L = [Z 1 / 9 ]. 

To obtain a lower bound on tr(i?^ z d) we first use the lower bound Jl3], |ll) 

on the correlations, namely — const[iV J p^/ 3 ] 1 ^ 2 , to reduce it to a radial problem. 
Using the fact that Z/r > (p^ F (r) for r > it follows from this that 

tv{H^ z d) 

L — X oo 

>^2q{2l + l)Y / e^ l ^ z + Y / q{2l + l)J2en,i,e,z-D{p TF ,p TF )-constzi. 

1=0 n l=L n (39) 
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Note that (39) arises from a relatively simple lower bound calculation. Part of 
the proof of the Scott conjecture amounts to proving that the right hand of (39) 
is accurate to o(Z 2 ). This proof was carried out in |23| (see also || 24|). We 
are not rederiving the Scott correction for the energy, and it is not necessary for 
us to do so here. 
Define 

' \ n > 
n < ' 



9(n) 



Since the eigenvalues of the perturbed problem (e ^ 0) are equal to the unper- 
turbed one (e = 0) except for I = Iq, we get the inequality 



lim sup X(Z) 



< lim inf lim sup 

«V Z-»oo 



llo + 1) W ' == l °' e ' z) 0(L - l ) 
eZ z 

tr(ff; , ,z)- - tx(Hi ^,z) 



-9(l - L) + constZ 



(40) 



eZ 2 

Because L eventually becomes larger than the fixed Iq, 

limsupA(Z) < q(2l + 1 ) lim inf lim sup °'° '— °" 

tr(i7, H n ,)_ - tr(H, H f 
= g(2f + l)limiaf °' oa V l °' eA ' , (41) 

e\0 e 

where the last equation holds since because of the scaling of hp z and hi 0:<L: z- 
Therefore 



lim sup X(Z) 

Z^oo 



< g(2/ + l)liminfy e ^"'°- 1 (42) 



e\0 



= <?(2Z + 1) Vliminf^iM r^Mil (43) 



,2 „i? 



(44) 



To exchange the limit e \ with the summation in ( (42| ) we use Lemma ^], 
which provides a summable majorant for the series that is uniform in e and 
thus allows us to fulfill Weierstrafi' criterion for uniform convergence. Finally, 
to deduce (|44|) from ( |43] ) we use the fact that the one-dimensional delta potential 
is a relatively form bounded perturbation, i.e., defines an analytic family in the 
sense of Kato. 

To obtain a lower bound for \{Z) we pick e negative instead of positive and 
take the limit limsup c ^ an d lim inf z^oo instead of lim inf ^0 and limsup^^^. 
Repeating the same steps gives the same result except for reversing the inequal- 
ities, thereby yielding the same bound (^J) from below. This establishes the 
first claim of the theorem. 
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2. Proof of the weak convergence: Because of the linear dependence of the 
right and left hand side of (|l2|), it suffices to prove the claim for the positive 
and negative parts of V separately. Thus we may - and shall - assume that 
V is positive. We can now roll the proof back to the previous case as follows: 
First we pick Z large enough so that Iq < L. It is convenient now, to replace e 
be e/a in order that the right side of j32|) in Lemma § is uniformly bounded in 
a and e for all a £ (0, oo) and for |e| < n/8. Then we integrate the inequality 

tr(ff, ff n 7 )_ -tr(H" , „)_ 

a 2 pi (a) < \ 7 2 6{L-l) + const/{(e/a)Z™), 

{e/a)Z (45) 

i.e., 

ap; (a) < ^^-^ 0(L - I) + coast /(eZ™) 

eZ (46) 

against V(a) from to oo. Thanks to Lemma [j] the right of side of ([46]) is 
bounded by const a and hence the integral is finite. Next we write out the 
traces appearing in (|44|) in terms of the eigenvalues and then use Lemma || to 
provide a bound that is summable (over the eigenvalues) and integrable (from 
to oo), if |e| < 7r/8. This bound is uniformly bounded in e, and so, by dominated 
convergence, we can take the limit e \ term by term. Using the result (|TT|), 
which we established above, Equation (|l2| ) is now verified. | 

Proof of Theorem [|: As was the case in the proof of Theorem we shall 
assume that W and V are nonnegative. For Part 1 we proceed as for Theorem 
H and define z as in (|34|), but with n(Z ) replaced by W(lo). First we treat 
the case W(u>) — 1. We follow the proof of Theorem |l| up to Equations ( |36"| ) 
and @ (with <5 Mo replace by W). The n we obtain analogously 

limsupA(Z) (47) 



Z^c 



< lim inf lim sup 

e\0 z^oo 



<l(2l + 1) 1 ' eZ2 " 9{L - 1) (48) 
.1=0 6 



Erm , -,^ T (Hi,o,z)- - tr(#/. e .z) -o-i- 
g(2Z + 1) — ^ — '-^ — -3— - — 9(1- L) + const Z «e 



(49) 



Vg(2Z + l)liminflimsup— '' £ ' Z (50) 

V g(2Z + l liminf — — (51) 

z=o 



a 2 £ lim. inf (a) (52) 



z=o 

,2 „H 



= a 2 p H (a). (53) 

To obtain ( |49| ) we use inequalities (38) and (39). To obtain ( p0[ ) we use the 
fact that Lemma [l] provides a majorant uniform in e and Z which is absolutely 
summable with respect to J2l=o ( l(^ + -0> i- e -' fulfills the Weierstrafi criterion 
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for uniform convergence (or the hypothesis of Lebesgue's dominated convergence 
theorem), and therefore allows the interchange of the limit and the I summation, 
and that the second sum tend term by term to zero. To obtain ( pl| ) we use 
the fact that the eigenvalues of the bare problem scale like Z 2 . Finally, the 
convergence result of Theorem 2 was used to obtain (|5|) . 

To obtain a lower bound we pick e negative instead of positive and take the 
limit limsup e ™ and liminfz^oo instead of liminf e ^o and limsup^^^. Re- 
peating the same steps gives the same result except for reversing the inequalities 
thereby yielding the bound from below. 

Let W now be a general bounded, measurable function on the unit sphere 
which we may - according to the remarks in the beginning - assume to be 
positive. We take ||W||oo=l- 

Let us try to imitate the steps ([!?]) to ( |53| ) . As before we are faced with esti- 
mating the eigenvalues of the one-body operators Hf z := —A — Z/\.\ — eZWS^ 
and Hj z = —A — ip^ F — cZW5<l but unlike the previous case these cannot be 
simply indexed by the angular momentum I when e ^ 0; indeed the one-body 
operators cannot be reduced to a direct sum of radial Schrodinger operators as 
in ([36]) and (|37]). However, the eigenvalues are real analytic functions of e and 
we can label the eigenvalues by the lvalue they have when e tends to zero. In 
short, the only change needed in (|47| ) to ( |53| ) is to replace (21 + l)e^ e z by the 
sum of the eigenvalues in the multiplet of H^ z that converge to z as e 
tends to zero. Since W is bounded by 1, all our previous bounds for eigenvalue 
differences (Lemmata |l| and ^) continue to hold and we are finally led to the 
lim inf in (^) . 

The crucial point is this: Even if W is not spherical symmetric, the sum of 
the eigenvalues in any multiplet is rotationally invariant to first order in e in 
the following sense. The only property of W that matters - to first order - is 
the average W averag e ■= (47r) _1 / W(u))du. 

Reversing the sign of e again gives the lower bound. 

2. Proof of the weak convergence: The proof can be rolled back to the 
previous case as follows: First we assume that V is spherically symmetric and 
integrate the inequality 

a 2 p z (auj) 



< 



<7(2Z + 1) l '°' Z 2 V l ' e ' Z> 9{L - I) 

.1=0 e 

+ f) q(2l + 1) tr(g|.o,z)--tr(g,,,,z)- g(f _ L) 



1=0 



const/ (eZ 24 ) 



against aV(a) from to oo. Observe that because of Lemma ^ the summand of 
the sum on the right of side of this integrated inequality is uniformly bounded 
by 

9 f°° 

(i + im + i)J v{a)ada 

which, when multiplied by (21 + 1), is summable. Again, the same argument 
holds when expressing the traces as sum over eigenvalues. Thus we are allowed 
to take the limits term by term for the differences of the eigenvalues giving the 
desired result as above. 

The extension to the non-spherical case is as in Part 1. | 
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5 Extensions to Molecules 



The ground state energy of a neutral molecule with nuclear charges Z\ = 
Xzx, Zk = Azk and positions of the nuclei at ...,!lRr- is given as 

E(N, Z) = inf{inf a(H N ^\R G R 3K } (54) 

where 

N ( K Z K \ N 1 Z K Z K , 

h nXR = E - E + E ^— ^ + pt=% 5) 

self-adjointly realized in i]Ar. Here Z denotes the AT-tuple [Z\, Zk) and i? 
the 3-FT-tuple (■EH 1 , ...,3tjf). We also set z := (zi,...,z^). Solovej S showed 
recently that for arbitrary but fixed z and N = Z\ + ... + 

A' 

£(7V,Z) = ^£(Z K ,Z K ) + o(Ai) (56) 

re=l 

holds as A tends to infinity and that the minimizing inter-nuclear distances are 
of order A~ 5 / 21 or bigger. These results imply among other things not only that 
the atomic Scott correction and Schwinger correction implies the molecular one 
but allows us to generalize Theorem ^ as well: The molecular density in the 
vicinity of each nucleus converges in the sense of Theorem 1 to the hydrogen 
density at each of the centers. Our precise result is: 

Theorem 3 Assume that E(N, Z) as defined in $5^ ) is equal to 

w£{va£a(H N g R )\R G M 3K , Vi< k<k <<aH«k - K\ > R ■= constA 7 } 

(57) 

with 7 > —1/4. Assume N = Z\ + ... + Z^, Z\ = \z\,...,Zk — Xzk with 
given fixed z\,...,zk- Furthermore fix kq G 1,...,K and pick a sequence of 
ground state density matrices d c of H N g ^ with densities pd c ■ Define p\, K „ (r) := 
p\((v- ?r\ Ko )/X)/\ 3 . Finally assume We L 2 (S> 2 ). Then 

W(u)px,*(™) -> qp H (r) [ W. (58) 



s 2 

Proof. First note that by suitable relabeling we can always assume that kq = 1. 
Set H € N ^ £ '■— H N g — Ylv=i e ^( r — ^i)- Because of ( |j6"| ) it suffices that 

N 

tr(dW N - R ) > tr(if! - eU x )+ J>(ff K )_ - D(p%,p%) ~ constZ 2 ^ 

for some positive S and an approximate ground state d. To this end let us 
introduce the localizing functions 

v K (x) :=cos(iP(\t-M K \/R)) 

where tp(t) is some continuous, piecewise differentiable, monotone decreasing 
function which vanishes, if t < ^, and which is n/2, if t > i. Note that the 
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supports of these functions have at most finitely many points in common because 
R is the minimal nuclear distance. We also define 



K 



v :- 



\ 



Now pick the density p(v) := Y^, K =i Pz F (\ v ~^-n\) an d denote the one-particle 
density matrix belonging to d by d\. - Note that tr d\ = N . - By the correlation 
inequality and the localization formula using the above decomposition 

of unity we have 



> tr < 



K 

k' = 1 
A" 



E 



z K z K , 



- p K - eU z 



— const A 3 



K 



^ tr { E u « 



A" 



-A3-^^(|.-^5|)-eC/ z 



AT 



A 



^Igradu^HooiV- ^ p K ) - constA^ . (59) 



t=0 



' = 1 



In (59) we used the spherical symmetry of <pi, ipx to show that D(p K ,p K i) < 

Now pick any arbitrary pair of different indices k, k' G {1, ...,K}. On the 
support of of v K we have 



<PVj\*-n*>\)< 



2 2 3 4 ^ 2 
q 2 (R/2) 4 



where we use the fact that the Sommerfeld solution of the Thomas-Fermi equa- 
tion is a pointwise upper bound of the Thomas-Fermi potential ( |1^| , Section 
V.2). Thus on the support of v K we have 



K 



i? 3 ' 



For the derivative of the v K , which governs the localization error, we have 
the following uniform estimate: w^^L | gradu re | 2 = |^'(|t — 5H K |)| 2 /i? 2 = 4n 2 
where, for definiteness, we picked if) to be the linear functions interpolating 
between and ^ on the interval [4, i]. Note that outside the annuli of thickness 
R/2 centered at the nuclei the derivatives vanish, in fact, whereas in these annuli 
the bound is actually an equality. 



15 



> 

K=l 



This yields 

^K.Z^ 

* f \f At: 2 2 4 3 4 ^ 2 \ , 1 „. J 

K— 1 ^ ^ 

> inf <r(fT| liZl ) + £inf ffOffz,,.^ - const (A + A 5 / 3 ) 

if 

> inf +^inf *(H Zk , z J - o(X 2 ). (60) 



At — 2 



Combining this with Solovej's upper bound reduces the converge question to 
that of the one-center casc.l 



A Appendix: Facts about the Atomic Ground 
State Energy 

According to |22| we have 

Ez.z < E TF (Z, Z) + § Z 2 + constZ^, (61) 
8 

and according to |^| (see also Q and Hughes ||) 

L-l 

Ez.z > Y,q(2l + l)tr(H l H 0>z )_ 
1=0 
00 

q(2l + 1) tr {H l>0 ,z)_ ~ D{ Ptf ,Ptf) - constat 



l=L 



> E TF (Z, Z) + \z 2 - constZ^r l og Z (62) 
8 

with L = [Z'b]. Combining (^lj) and (^32|) gives 

L-l 

Sz.z = 5>(2J+l)tr(fff 0!Z )_ 

00 

E47 
Q(2Z + 1) tr (fr,,o,^)_ - D( PTF , PTF ) + 0(Z*). (63) 



l=L 
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